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Abstract. - We introduce a bond percolation procedure on a D-dimensional lattice where two 
neighbouring sites are connected by N channels, each operated by valves at both ends. Out 
of a total of N, randomly chosen n valves are open at every site. A bond is said to connect 
two sites if there is at least one channel between them, which has open valves at both ends. We 
show analytically that in all spatial dimensions, this system undergoes a discontinuous percolation 
transition in the N — > oo limit when 7 = j^jl crosses a threshold. It must be emphasized that, in 
contrast to the ordinary percolation models, here the transition occurs even in one dimensional 
systems, albeit discontinuously. We also show that a special kind of discontinuous percolation 
occurs only in one dimension when N depends on the system size. 
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Percolation transition is one of the most studied [T] crit- 
ical phenomena in non-equilibrium statistical physics. It 
is usually modeled on a lattice [5] where bonds connecting 

] the neighbouring sites form independently and randomly 
with probability p. In two and higher dimensions, a con- 
tinuous phase transition to a percolating state having an 
infinitely large connected cluster occurs when p crosses a 
threshold value p c > [TJ|3]. One dimensional systems, 

■ on the other hand, do not show percolation transition as 
these systems cannot have a percolating configuration for 
any p < 1. The trivial fixed point p c — 1, which can be 
approached only from the non-percolating region, shows 
signature of a continuous phase transition. In this article 
we propose a generic bond percolation procedure which 
shows a percolation transition, even in one dimension. It 
turns out that the transition is discontinuous. 

Discontinuous percolation transition has been a subject 
of great interest since the study of Achlioptas et. al. [4]. 
They proposed a non-local bond percolation procedure 
which initially inhibits formation of a single large cluster. 
When applied to a fully connected graph, it apparently 
results in a discontinuous transition, which was named as 
explosive percolation transition (EPT) . Different aspects of 
such explosive transition resulting from Achlioptas growth 
process on several other graphs have been studied in large 
number of following works 016]. In recent studies [7] it 
has been shown that EPT can also be obtained from per- 
colation rules other than Achlioptas process. In fact, it 
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Fig. 1: (Color online) The ID model with N = 4 channels : 
each site has N = 4 valves (discs), of which n — 2 are open 
(blue) and the rest are closed(red). Existence of at least one 
channel which has open valves at both ends (shaded), repre- 
sents a bond between the corresponding sites. 



has been argued that certain specific rules [5] of adding 
new bonds, which depend on the properties of the existing 
clusters, are sufficient for a system to show EPT. However, 
discontinuous nature of EPT has been questioned recently; 
this originated a long debate [9] regarding the order of the 
transition. It has been argued that explosive percolation 
is, in fact, a continuous transition [10 with a small crit- 
ical exponent /3 which only appears to be discontinuous 
in numerical studies. A percolation process is proposed 
in the following which unambiguously shows a discontin- 
uous transition. Both the critical point and the jump in 
the order parameter there in the thermodynamic limit are 
calculated exactly. 

In this Letter we introduce a local percolation proce- 
dure on a D-dimensional lattice where every pair of neigh- 
bouring sites has N different channels joining them. Each 
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channel contains a valve at every site which can either be 
open or closed independent of the other valves. Two sites 
are connected by a bond if there exists at least one channel 
which has open valves at both the sites. Clearly each bond 
appears randomly, independent of other bonds and clus- 
ters. This system shows a discontinuous transition in the 
thermodynamic limit when the relevant tuning parameter, 
namely the number of open valves n at each site, crosses 
a threshold value. The model could be solved exactly to 
locate the transition point. We show that, the transition 
is discontinuous in all dimensions, including ID. 

First let us describe the model in one dimension; its 
extension to higher dimensions is straight forward. The 
sites labelled by i = 1, . . . , (L + 1), are connected to their 
neighbours by A channels (see Fig. [IJ. Each channel 
contains a valve at every lattice site i that can either be 
open or closed. Out of a total of N valves at every site, 
n are chosen randomly and opened. The neighbouring 
sites of the lattice are said to be connected by a bond 
if there exists at least one channel between them which 
has open valves at both ends. The model can be recast 
into a simpler form by associating a set Si of n integers, 
randomly chosen from a larger set {1,2,... TV}, to each 
lattice site i. In this picture, a bond is said to connect two 
neighbouring sites i and j if Si (1 Sj is not null. 

This system, being one dimensional, is percolating only 
when all the L bonds are present. Obviously percolation 
is not possible when the number of open valves n at each 
site is zero. Again, the system is surely percolating for 
n > A/2, as in this regime every pair of sites has one or 
more common valves which are open. Our aim here is to 
find if the system is percolating for any non-zero n smaller 
than A/2. 

The principal quantity of interest is the percolation 
probability Pl, which is defined as the probability that 
an arbitrary lattice site belongs to the spanning cluster 
[US]. In other words, Pl is the average number of sites 
belonging to the spanning cluster and plays the role of 
the order parameter as it is nonzero only in the percolat- 
ing regime for a thermodynamically large system. For the 
usual bond percolation, in two and higher dimensions, Pl 
vanishes continuously at the critical point for D > 2. It 
is worth mentioning that Pl should not be confused with 
TIl, the probability that there exists a spanning cluster in 
the system. In fact, in the thermodynamic limit IIoo for 
all spatial dimensions jumps from to 1 as the connection 
probability p is increased beyond p c . The fact that IIoo is 
discontinuous, is sometimes used p] to locate the exact 
transition point. 

Now let us calculate TIl for this multi-channel model in 
one dimension. 
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where q is the probability that two neighbouring sites are 
not connected by a bond. TIl vanishes in the thermody- 
namic limit L — > co for any q > 0, which corresponds to 
the connection probability 1 — q = p < 1 . 




Fig. 2: (Color online) (a) Phase diagram in the ot-u plane: 
u c — g^ 1 (a) is the critical line. For a > In 2 transition occurs 
trivially at v c = ~. (b) Analytical expression of Pl as a func- 
tion of v (lines) for a = 0.01 and N = 400, 800 and 1600 (cor- 
respondingly L = e aN ) are compared with the numerical sim- 
ulations (symbols). The discontinuity, occurring at v c — 0.095 
is more prominent for larger N. 



For this model q can be calculated as follows. Since n 
valves can be chosen from the total N in possible ways, 
the probability that k open valves are common between 
any two neighbouring sites is 



Q k = 



r^N—nf-in 



(2) 



Thus, a bond between any two neighbouring sites i and 
i + 1 is absent with probability q — Qo, which can be 
expressed as a function of A and v = (the density of 
open valves at each site). Using Stirling's approximation, 



(1 



exp [-g{v)N] 



with g{v) = — In 
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Note that the probability p = 1 — q of connecting neigh- 
bouring sites by a bond, here, is independent of other 
bonds and clusters. 

From Eq. © it is evident that q is finite, though ex- 
ponentially small, for any given N. Thus, like other one 
dimensional models, here too, one cannot have a perco- 
lating state in the thermodynamic limit L — > oo for any 
arbitrary value of n < N/2. However, for any given L, if 
one lets N — > oo, then q vanishes, resulting in a percolating 
state for all n > 0. To explore the possibility of transition 
at a non-trivial n, let us couple A to the system size L. 
Since q is an exponentially decaying function of N, it is 
suggestive that one takes A = O(lnL). Let the thermo- 
dynamic limit be taken in such a way that L approaches 
oo along with A, whereas their ratio a = =th& remains 
fixed. From Eqs. (P and ©, now, U L = (1 - cL~ 9 MA*) £ 



(1) where c 



lim TIl 



in the thermodynamic limit, 

1 for g(v) > a 
e~ c for g{y) = a 
for g{v) < a. 



(5) 
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Fig. 3: (Color online) (a) The order parameter Pl as a function 
of v for N = 400, 800, 1600. Pl approches a step function at 
v c = 0.095 as N — > oo. (b) The jump in the order parameter, 
A (Eq. pp. for different values of 8v approches 1 as L is 
increased. In both the figures, the system size is L = e aN with 
a = 0.01. 

Thus, transition from a non-percolating state (IIoo = 0) 
to a percolating state (IIoo = 1) occurs at v c — g~ 1 (a). 
The corresponding phase diagram is shown in Fig. EJa). 
Note, that for a > In 2, the transition occurs at the trivial 
value v c = 1/2, which corresponds to half of the valves 
(n = N/2) being open at every site. In Fig. HJb) we have 
shown IL L versus v for three different values of N with 
a = 0.01. Symbols therein represent the same obtained 
from Monte-carlo simulations. Clearly LI^ approaches the 
step function Q(v — v c ) as L — > oo. 

In one dimension, the existence of a spanning cluster 
(i.e., when IIoo = 1) implies that all the sites of the lattice 
belong to that cluster; which in turn implies that in the 
percolating regime, the order parameter is P^ — 1. This 
clearly indicates that the percolation transition is discon- 
tinuous. In Fig. Ela), we have shown P L for L — e aN with 
a = 0.01 and N = 400,800, 1600. To show that P L , like 
ILl, also approaches Q(u— z/ c ) as L — > oo, we estimate the 
jump in the order parameter P L across the critical point 
v c as 

A(L; 5v) = P L (y c + Sv) - P L {v c - 5v). (6) 

In Fig. [31(b) we plot A as a function of L for 8v = 
0.0075,0.0125,0.0225. The monotonic increase of A with 
L is a clear indication that has a discontinuity at v c . 
For any 5v > we have lim A(L; 5v) = 1. It is evident 

from the figure that this limit is approached extremely 
slowly. 

Until now we have been discussing possibility of bond 
percolation transition in this one dimensional system 
where both N and L approach infinity keeping a — 
fixed. If instead n is related to N, one can obtain another 
percolation transition. Note, that q in Eq. (j3]) vanishes 
exponentially in the limit N — > oo, only when g(y) > 0. 
However, if g(v) — > slower than 1/N, q can approach 1 
even in the TV — > oo limit. Let us discuss this scenario in 
details. 

It is clear from Eq. (0} that g(y) vanishes only at v = 0. 
To the leading order in v we have g(y) = v 2 and c = 1, 
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Fig. 4: (Color online) $(7, N) versus 7 (main figure) and Pl 
versus 7 (inset) for N = 10 2 , 10 4 , and 10 12 . Simulation results 
for 4>(7,iV) with L = 4 (circles) and L — 9 (crosses) are also 
shown in the main figure. The inset shows Pl for L — 100 as 
a function of 7 for different N. 

thus g{v)N = n 2 /N. Let us take n = N~< with < 7 < 1, 
so that v approaches zero in the large N limit. This results 
in 

q = l-p = cxp(-N 2 "<- 1 ). (7) 

Now, in the N — ¥ 00 limit, the connection probability 
p = 1 — q is 

r for 7 < 1/2 

lim p = I 1 - e" 1 for 7 = 1/2 (8) 

N-^-oo 

[l for 7 > 1/2. 

Thus, again, as 7 crosses the threshold value j c = |, i.e. 
when the number of open valves at each site n is increased 
beyond y/~N, both P L and ITl in the thermodynamic limit 
jump discontinuously from to 1, resulting in a discontin- 
uous percolation transition. 

It is difficult to simulate this model with large L for 
the reason that the limiting values of p [Eq. ©] are ap- 
proached extremely slowly. Accordingly, one needs un- 
reasonably large value of N to see the transition at j c . 
For example one needs N ~ O(10 30 ) for a system of size 
L = 100. This difficulty can be avoided if we measure 

$(7,iV) = i- 1 lnn i =ln(l-g) = In (l - e^ 2 "^ 1 ) , (9) 

which is independent of L. A distinct signature of this 
discontinuous transition is that the curves $(7, N) versus 
7 for different values of N intersect at 7 = 7 C as at this 
point $(7c, N) = ln(l — e _1 ) is a constant. In Fig. |4]we 
have plotted $(7, JV), obtained numerically [symbols] for 
systems of size L — 4, 9 , as a function of 7 with two differ- 
ent N = 10 2 , 10 4 and compared those with Eq. © [solid 
lines] . An excellent match between these two for different 
is suggests that the transition is present for systems of 
all sizes. It is only that the required numerical accuracy 
for large systems is hard to achieve in affordable computa- 
tional time. A plot of Eq. © for N = 10 12 , which could 
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Fig. 5: (Color online) (a) The percolation probability Pl of 
usual bond percolation on a 100 x 100 square lattice is plotted 
against p, along with 7 versus p [calculated from Eq. 0] for 
our model in 2D with N = 2048. For any value of 7, which 
corresponds to a unique p (thus Pl) one can read out Pl using 
(0 [following the arrow]. The resulting Pl for N = 2048 and 
10 12 are shown as solid lines in (b). Symbols therein correspond 
to Pl obtained from direct numerical simulation of our model 
on a 100 x 100 lattice with N = 2048. 

not be supplemented with the corresponding numerical re- 
sults, is shown in the same figure to demonstrate that the 
transition, in fact, occurs at 7 C = i. 

Next, let us look at the order parameter Pl- Since 
p(7c) = 1 — 1/e for large N, one expects that the plots 
of IIl versus 7 for different values of N (reasonably large) 
would intersect at j c . This feature is clearly visible in the 
inset of Fig. 01 where we have plotted Pl for a system of 
length L = 100 as a function of 7 for N = 10 10 , 10 20 , 10 30 , 
and 10 60 . For these large values of N we have calculated 
Pl, the average size of the largest cluster, by connecting 
the sites with probability p obtained directly from Eq. (J7J . 
It is evident that the order parameter is discontinuous in 
the N — > 00 limit. 

Let us summarize the results obtained till now. We 
show that a discontinuous transition can be obtained in 
this one dimensional system as N and L go to infinity, 
in two different ways. (A) Both a = and v — 4r 
are tuned to obtain the transition, yielding a non-trivial 
phase diagram (Fig. Ufa)) in the a-v plane. (B) 7 = 
plays the role of the control parameter resulting in a phase 
transition at j c — ~. The natural question to ask next is, if 
any of these discontinuous transitions is present in higher 
dimensions. 

The model can be extended to two dimensions in a 
straightforward manner; we choose to work on a square 
lattice. The neighbouring sites of this L x L lattice are 
joined by N channels. Correspondingly, there are N valves 
at each site which can either be open or closed. An open 
valve at a site allows possibility of connection to its neigh- 
bours in both vertical and horizontal directions. As before, 
two neighbouring sites of the square lattice are connected 
by a bond, only if there exists at least one channel which 
has open valves at both these sites. Thus, the probability 
q that two neighbouring sites are not connected by a bond, 
when n randomly chosen valves are opened at each site, is 
again given by Eq. ([3]). 



In the usual bond percolation problem [I] , as mentioned 
earlier, the percolation transition is governed by an or- 
der parameter Pl which is defined as the probability that 
a randomly selected lattice site belongs to the spanning 
cluster. The probability that a spanning cluster exists, 
i.e., IIl, changes discontinuously across p c in the L —¥ 00 
limit, as in case of ID. However, unlike one dimensional 
systems, the exact formula for 11^ is not known in higher 
dimensions. Henceforth we concentrate only on the order 
parameter, denoted by Pl{p) for a L x L square lattice. 
In the thermodynamic limit L —± 00, Pl{p) vanishes con- 
tinuously at the critical threshold p c = 1 — q c = i . 

It is well known that the percolation probability Pl{p) 
of a finite two dimensional system differs from (p) only 
by a correction factor which is negligibly small for large 
L. This indicates that a discontinuous transition similar 
to case (A) cannot be obtained in two dimension just by 
rendering N a function of L. On the other hand, for any 
arbitrary N, the percolation probability Pl for this model 
on a L x L lattice can be obtained from Pl as 

PL=P L (p=l-q), (10) 

where q is given by ([3]). So, the usual bond percolation 
that occurs as a continuous phase transition at p c = 1/2 = 
q c is also expected here for any given N when v increased 
beyond a critical value v c which is a solution of e 9 ^ c ^ N = 
2c(y c ). This continuous phase transition is similar to the 
usual bond percolation transition [3] on a square lattice. 

Now we turn our attention to case (B) where 7 = 
is used as a tuning parameter. Here, the connection prob- 
ability p has two distinct limiting values [see Eq. ©] 
as N — ¥ 00. So, when 7 is varied continuously in the 
range [0, 1] the connection probability p jumps from to 
1 at 7 C = ^, resulting in a discontinuity in the percola- 
tion probability Pl across j c . Thus, this model shows a 
discontinuous transition also in two dimensions. 

The order parameter Pl can be calculated from Pl{p) 
using Eq. (fTOf . As the analytical form of Pl{p) is not 
known, we first obtain the same numerically by simulating 
the usual bond percolation problem on a L x L square 
lattice for < p < 1. Pl can be read out from this data 
using Eq. (I10[) where p = 1 — q is found from Eq. (J7J for 
any given 7 and N . This procedure is illustrated in Fig. 
(Ha) for a 100 x 100 lattice and N = 2048. The resulting P L 
is shown in Fig. OUb) as a solid line. The symbols therein 
correspond to the same obtained from direct simulation of 
a 100 x 100 square lattice having N = 2048 channels. It is 
clear from this figure that the transition occurs explosively 
near j c = 1/2. Pl for N = 10 12 , obtained using the above 
procedure, is also shown in Fig. E^b) to illustrate that, as 
expected, the transition point shifts towards j c — 1/2 as 
N -> 00. 

In contrast to the usual continuous bond percolation 
transition, the spanning cluster appears suddenly in case 
of a discontinuous transition, as one approaches the tran- 
sition point. In Fig. [6] we have shown the snap shots of 
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Fig. 6: (Color online) Typical configurations of this model (top 
panel) for N = 2048 at 7 = 7 C — e is compared with that of 
the usual bond percolation (bottom panel) at p = p c — e. In 
both the cases 50 large clusters on a 500 x 500 square lattice 
are shown. 

the typical configurations (50 large clusters are shown) of 
our model on a 500 x 500 lattice near the transition point 
taking 7 = 7 C — e and compared the same with the usual 
bond percolation at p = p c — e. Clearly, when e — > 0, 
the size of the largest cluster (shown as blue) grows faster 
in the former case, indicating the explosive nature of the 
transition. 

It is straight forward to study these models on other 
kinds of lattices in two or higher dimensions. The explo- 
sive transition that occurs when the number of open chan- 
nels is varied as n = N* 1 , is quite generic. In fact the limit- 
ing values of the connection probability p [Eq.©] changes 
from being to 1, at 7 C = h when number of channels 
N — > 00. Corresponding percolation probabilities change 
from to 1 at 7 C resulting in an explosive percolation in all 
dimensions. As discussed earlier, the other scenario where 
an explosive transition was obtained by tuning v = jr and 
a = tt, is specific to only one dimensional systems. 

In summary, the procedure we introduce in this arti- 
cle, though a simple A^-channel extension of the usual 
bond percolation, shows discontinuous percolation tran- 
sition in all dimensions. Unlike the usual bond perco- 
lation transition, which is not possible in ID and occurs 
as a continuous phase transition in higher dimensions, this 
model shows discontinuous transition in all /^-dimensions, 
including ID. The neighbours here are connected by N 
channels, each having an operating valve at either ends. 
Of these N valves at each site, n are open. The neigh- 
bouring sites are said to have a connecting bond if one or 
more channels joining them are open at both ends. We 
show that, when 7 = crosses a threshold value 7 = 5 
this system percolates abruptly in all spatial dimensions. 
It is rather surprising that a percolation transition occurs 
even in ID at a non-trivial value of the tuning parameter. 
The reason lies with the fact that in the N —> 00 limit, 
the connection probability p changes discontinuously from 



to 1, even though the tuning parameter 7 is varied con- 
tinuously. 

Along with the above transition, another explosive per- 
colation transition occurs especially in one dimension, 
driven by two parameters a — and v = j,. Occur- 
rence of this transition owes to the fact that the percola- 
tion probability Pl in ID explicitly depends on the sys- 
tem size L, unlike in higher dimensions where system size 
merely appears as a correction term. In fact a similar tran- 
sition is possible in the usual one dimensional bond perco- 
lation scenario if one uses a L-dependent q = L~ z , i.e. the 
connection probability p = 1 — L~ Z . Clearly, in the thermo- 
dynamic limit, the percolation probability Pl = (1 — q) L 
jumps from to 1 at z c = 1. 

We conclude with a few comments on the differences 
of this model with the explosive percolation. The discon- 
tinuous change in the size of spanning cluster have been 
reported earlier by Achlioptas et. al. [4] and several fol- 
lowing works 5 - 7 under the name of explosive percolation 
transitions. Later studies, however, indicate that these 
transitions are in fact continuous |10) . with an unusually 
small critical exponent. The percolation transitions re- 
ported here are analytically proven to be discontinuous. 
In all the studies of explosive percolation, the connection 
probability is allowed to evolve with the clusters, which 
facilitates the sudden formation of the spanning cluster 
resulting in an abrupt change in order parameter. In con- 
trast, the connection probability of the model studied here 
does not at all depend on the existing clusters; every bond 
appears with the same probability. The resulting discon- 
tinuous percolation is truly an emerging behaviour. 
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